UNIFORM ASYMPTOTICS FOR DISCRETE ORTHOGONAL 
POLYNOMIALS WITH RESPECT TO VARYING EXPONENTIAL 
WEIGHTS ON A REGULAR INFINITE LATTICE 



PAVEL BLEHER AND KARL LIECHTY 

Abstract. We consider the large N asymptotics of a system of discrete orthogonal poly- 
nomials on an infinite regular lattice of mesh i, with weight e ~ NV ( x ) ; where V(x) is a real 
analytic function with sufficient growth at infinity. The proof is based on formulation of 
an interpolation problem for discrete orthogonal polynomials, which can be converted to a 
Riemann-Hilbert problem, and steepest descent analysis of this Riemann-Hilbert problem. 

in 

de 1. Introduction 

■i— > 

a 



For a given N G N, introduce the regular infinite lattice, 



j n = s Xk,N = , k e z| . (1.1) 

We consider polynomials orthogonal on with respect to the varying exponential weight 

oo 

where V(x) is a real analytic function such that, for some e > 0, V has analytic extension 
into the strip 

\lmz\<e, (1.3) 

<^> and satisfies the growth condition 

Re V(z) 



, . — 7- +oo as Id oo, llni^l < e. (1.4) 

X' iog(kr + i) 

More specifically, we introduce the system of monic orthogonal polynomials, 
P n (x) = x n + pn^.-ix™" 1 + . . . + Pno, n = 0, 1, . . . , 

such that 

^2 P m {x)P n {x)w N {x) = h n S mn , (1.5) 

for some normalizing coefficients h n . Existence and uniqueness of this system of orthogonal 
polynomials is guaranteed by condition (1.4). These orthogonal polynomials satisfy the 
three-term recurrence relation 

xP n (x) = P n +l(x) + f3 n P n {x) + J 2 n P n -l(x). (1.6) 
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We will explore the asymptotics of the quantities 7„, /3 n , and h n for n = N, N — 1, as well 
as pointwise asymptotics of the polynomials Pn(x) as N — > oo. 
The present work has the three predecessors: 

(1) the work [6] of Deift, Kriecherbauer, McLaughlin, Venakides, and Zhou, in which the 
large N asymptotics has been obtained for orthogonal polynomials with respect to 
varying exponential weights on the real line, 

(2) the work [1] of Baik, Kriecherbauer, McLaughlin, and Miller, in which the large N 
asymptotics has been obtained for orthogonal polynomials with respect to varying 
exponential weights on a lattice in a finite interval, and 

(3) the work [3] of Bleher and Liechty, in which the large N asymptotics has been obtained 
for orthogonal polynomials with respect to the varying exponential weight wn(x) = 
e - N (\x\-Cx) on the infinite lattice L^. 

Also, a very important ingredient comes from the work [11] of Kuijlaars, in which analytic 
properties of equilibrium measures with constraints are established. 

The asymptotic analysis of the polynomials Pn(x) in this work will be based on the 
Interpolation Problem for discrete orthogonal polynomials, which is introduced in the work 
[4] of Borodin and Boyarchenko (see also [1], [2], [3]). The asymptotic analysis of Pjv(x) 
will consist of three steps. The first step will be a reduction of the Interpolation Problem 
to a Riemann-Hilbert problem on a contour on the complex plane, which we accomplish 
following the general approach introduced in the paper [12] of Miller and in the monograph 
[10] of Kamvissis, McLaughlin, and Miller. The second step will be an application of the 
nonlinear steepest descent method of Deift and Zhou [8] to the Riemann-Hilbert problem 
under consideration, and the third and final step will be a derivation of the asymptotic 
formulae both for the orthogonal polynomials Pn(x) and for the recurrence coefficients. To 
apply the nonlinear steepest descent method to the orthogonal polynomials Pn(x) we need 
to study the corresponding equilibrium measure. 

2. Equilibrium Measure 

The significance of the equilibrium measure is that, as we will see, it gives the limiting 
distribution of zeros of the polynomial Pn(x). By definition, the equilibrium measure is a 
solution to a variational problem. Namely, let us consider the following set of probability 
measures on R 1 : 

M = {0<u<a, u(R r ) = l}, (2.1) 
where a is the Lebesgue measure, and let us introduce the functional 

H{y) = J J log 1 du(x)du(y) + J V(x)du{x), v G M. (2.2) 

The equilibrium measure minimizes this functional over some set of measures. In the case 
of continuous orthogonal polynomials, we minimize over the set of probability measures on 
the real line. However, in the case of discrete orthogonal polynomials, we must introduce 
the upper constraint, v < a, in order to account for an interlacing property of the zeroes of 
orthogonal polynomials. 

It is a general fact, (see, e.g. [15]) that for any system of polynomials orthogonal on the 
real line with respect to a real weight, the nth polynomial has n real distinct distinct zeroes. 
Furthermore, the zeroes of a system of discrete orthogonal polynomials satisfy an interlacing 
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property with regard to the location of the nodes of the lattice Ln, so that no more than 
one zero may lie between any pair of adjacent nodes. It therefore follows that, if we denote 
by Hn the normalized counting measure on the zeroes of the iVth orthogonal polynomial in 
our system, 

^(a, b) < b — a + — for any — oo < a < 6 < oo, (2.3) 
so that (j, < a, where \i = liniAr^oo p^. With this constraint in mind, we define 

E = inf H(v). (2.4) 

v<EM 

It is possible to prove that there exists a unique minimizer z/ , so that 

E = H(v ), (2.5) 

see, e.g., the works of Saff and Totik [14], Dragnev and Saff [7] and Kuijlaars [11]. The 
minimizer is called the equilibrium measure. 

The equilibrium measure uq is uniquely determined by the Euler- Lagrange variational 
conditions: there exists a Lagrange multiplier I such that 

/( > I for x G supp z/ 
log \x - y\dis (y) - V(x) < ^ , c , , (2.6) 

61 y| yyJ \ </ for x G supp (er — z/ ), V ; 

see the works [5] of Deift and McLaughlin and [7]. In particular, 

for , 6 supp ,„n supp (,-,„). (2.7) 

The equilibrium measure is possesses a number of nice analytical properties, as shown by 
Kuijlaars in [11]. We will use these analytic properties, so let us discuss the results of [11]. 
First, observe that the constraint u < a implies the existence of the density, 

= (2 - 8) 

We can partition R into the three sets 





= {x G 1 


::2 J 


f log 


\x - y\du (y) 


- V{x) = 1} 




I + 


= {i£l 


■ :2 J 


f log 


\x - y\dv (y) 


- V(x) > 1} 


(2.9) 


r 


= {i£l 


::2, 


f log 


\x - y\du (y) 


- V(x) < I}. 





The structure of the equilibrium measure is well described in the following theorem of Kui- 
jlaars, obtained in [11]. 

Theorem 2.1. (Kuijlaars) For any real analytic potential V(x) satisfying (1-4), the following 
hold: 

(1) The density p(x) of the constrained equilibrium measure u (defined in (2.5)) is con- 
tinuous. 

(2) The sets I + and I~ are both finite unions of open intervals. 

(3) The density p is real analytic on the open set {x : < p(x) < 1}. 
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(4) The density p has the representation 



p (x) = -yqt(x) for xei°ur, (2.10) 

where q± is the positive part of a function q\ defined on 1° U I which is real analytic 
on the interior of 1° U I . The function q\ is negative on I , so that 

p(x) = for are/", (2.11) 

and it is nonnegative on 1°, so that 

p(x) = -^(x) for xel°. (2.12) 

7T 

(5) The density p has the representation 

p(x) = 1 - -\/q£(x) for xeI°Ul + , (2.13) 

where q 2 is the positive part of a function q 2 defined on 1° U / + which is real analytic 
on the interior of 1° U / + . The function q 2 is negative on I + , so that 

p(x) = l for xel + , (2.14) 

and it is nonnegative on I , so that 

p(x) = 1 - -y/q 2 (x) for x € 1°. (2.15) 

7T 

Remark: It follows from equations (2.12) and (2.15) that 

-aAiO*) = i - -Vq&j for xel°, (2.16) 

7T 7T 

hence q\ and q 2 uniquely determine each other. 

Notice that, according to point (2) of this theorem, the connected components of /° are 
either closed intervals or isolated points. Since u has compact support, we can write 

/° = LJh-,/y (2.17) 

i=i 

where 

a, < Bj for j = 1, . . . , q 

n f i i ( 2 - 18 ) 

Pj<a j+1 for j = l,...,q-l. 

Notice that the intervals (— oo, cxi) and (B g , oo) are components of The interval (Bj, (Xj+i) 
for 1 < j < q is a component of either / + or / _ . We therefore adopt the notation 



Ay = |i G {1, . . . , q - 1} : (Bj, a j+1 ) C / | 

As = jj e {1, ...,q-l} : (Bj,a J+1 ) c / + |. 

We will call an equilibrium measure uq regular if the following hold: 

(1) qi and q 2 are non- vanishing on the interior of 1°. 

(2) /° contains no isolated points, so that otj < Bj for all j — 1, . 



(2.19) 
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Figure 1. The graph of the density function for a hypothetical equilibrium 
measure with q — 5. Bands are denoted by bold segments, saturated regions 
by dashed segments, and voids by thin segments. 

(3) If j e A v , then q'^j) ^ and q[(a j+1 ) ^ 0. Also, q[{ai) ^ and q[%) ^ 0. 

(4) If j e A s , then q' 2 {(5 3 ) ^ and q' 2 (a j+1 ) ^ 0. 

For the remainder of this paper, we will assume that our equilibrium measure is regular. In 
this case, the sets J°, I + and I~ are each finite unions of intervals, so that 

-oo < «i < /3i < a 2 < (3 2 < ■ ■ ■ < a q < (3 q < oo, (2.20) 
and we classify these intervals as follows: 

Definition: A void is an open subinterval j G A v , or one of the intervals (— oo, cti), 

(f3 q , oo). The union of all voids is I~ . 

Definition: A saturated region is an open subinterval (f3j, Oj+i), j G A s . The union of all 
saturated regions is I + . 

Definition: A band is an open subinterval (otj, (3j), j = 1, . . . , q. The union of all bands is 
the interior of 1°. 

Observe that p(x) = on any void (f3j, «j+i), p{x) = 1 on any saturated interval (f3j, ojj+i), 
and < p(x) < 1 on any band (aj, (3j), see Figure 1. In addition, at the end-points of any 
band, p(x) has a square-root singularity. Namely, if aj is a common end-point of a band and 
a void then as x — > +0, 

p( aj + x) = Cy/x(l + 0(x)), C = |gi(«j)| 1/2 > 0. (2.21) 
and if ctj is a common end-point of a band and a saturated region then as x — > +0, 

p(a j + x) = l-Cyfe(l + 0(x)), C = \q' 2 ( aj )\ 1/2 > 0. (2.22) 
Similarly, if j3j is a common end-point of a band and a void then as x — > +0, 

ptf, - x) = CVS (1 + 0(x)), C = \q'i(h)\ 1/2 > 0- (2-23) 
and if (3j is a common end-point of a band and a saturated region then as x — > +0, 

p(p j -x) = l-CVx(l + 0(x)), C =\q' 2 {P j )\ 1/2 >0. (2.24) 

In the next section we introduce the (^-function, which will be our means of exploiting the 
equilibrium measure. 
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3. The ^-function 
Define the g-function on C \ (—00, j3 q ] as 

g(z) = \og(z - x)dv (x) (3.1) 

J ai 

where we take the principal branch for the logarithm. Also, introduce the numbers Qj for 
j = l,...,q- 1 as 

r/3 q 

2n / p(x)dx for j G A v 

r ^ 

2ir / p{x)dx + 2naj + i for j G A s . 

Ja j+1 

Properties of g(z): 

(1) g(z) is analytic in C \ (—00, f3 q }. 

(2) For large z, 

^ g ■ [0* x j 

g(z) = logz - -j , 9j= I — dv (x). (3.3) 

j=i z J a i 3 



SI, = < 



(3) 



m = 1 ^| (3.4) 



is the resolvent of the equilibrium measure. 

(4) From (2.9), we have that 

{= V(x) + 1 for x G 1° 
> V(x) + 1 for x G 1+ (3.5) 
<V(x) + l for x G I~, 

where g + and g_ refer to the limiting values from the upper and lower half-planes, 
respectively. 

(5) Equation (3.1) implies that the function 

G(x) = g+(x) -g-{x) (3.6) 
is pure imaginary for all real x, and 

G(x) = 2ni / p(s) ds. (3.7) 

J X 

Thus 

r iSlj for pj<x< a j+1 , and j G A v 

\ iilj — 2-nix for f3j < x < <X/ + i, and j G A s . 

From (3.5) and (3.7) we obtain that 

/A 

2g ± (x) = V(x) +l±2ni / p{s)ds for x G 1°. (3.9) 

J X 
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(6) Also, from (3.7), we get that G(x) is real analytic on the sets I + , I , and on the 
interior of 1°. We can therefore extend G into a complex neighborhood of any interval 
of analyticity for p, and the Cauchy-Riemann equations imply that 

dG(x + iy) 



dy 



2tt P (x) > 0. 



y=0 



Observe that from (3.5) we have that 

G(x) = 2g + (x) -V(x) -I = 



-[2g4x)-V(x)-l], xel°. 



(3.10) 



(3.11) 



4. Main Results 

In this section, we summarize the main results of the paper. In order to do so, we must 
first introduce some notations. Introduce the numbers flj^ for j — 0, . . . , q as 

fNQj for j e A v 
ir + ML for j e A s 
<■ 2,N for ,=0 f 4 ' 1 ' 
^ for j = q. 

Q N — (S1i ; at, • • • , ^ 9 -i,at)- (4.2) 
q 

R(z) = l[(z-a j )(z-fi j ) (4.3) 
i=i 

and let X be the two-sheeted Riemann surface of genus g = q — 1 associated with ^jR(z) 
with cuts on the intervals (aj, We fix the first sheet of X by the condition 



and the vector 
Let 



^/R(z)>0 for z > (3 q (4.4) 

on the first sheet. 

Introduce the following homology basis on X. For any j e {1, • • • , q — 1}, let Aj be a cycle 
enclosing the interval (j3j,aj + i) (passing through the intervals (aj,j3j) and (a^+i, Pj+i)), 
oriented clockwise, such that the piece of Aj which lies in the upper half-plane also lies on 
the first sheet of X, while the piece of Aj which lies in the lower half plane also lies on 
the second sheet of X. Also for any j e {1, • • • , q — 1}, let Bj be a cycle enclosing the 
interval (ai,j3j) (passing through the intervals (— oo, «i) and (f3j, <x,- + i)), oriented clockwise, 

and lying entirely on the first sheet of X. Then the cycles (A 1: . . . , A q _ 1: B 1: -B ? -i) form 

a canonical homology basis for X. 

Now consider the the g-dimensional complex linear space Q of holomorphic one-forms on 
X, 

q = L = v^££}, (4.5) 

and the basis 

u = (wi, . . . ,w 9 -i) (4.6) 
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normalized such that 



/ u k = 5 jk . (4.7) 

J A; 



Notice that the basis oo is real. That is, for the basis elements 



CjkZ k 1 dz 



the coefficients Cjk are real. 

Now define the associated matrix of -B-periods as 



(T jk ), T jk = u k , j,k = l,...,q-l. (4.9) 



Since y/R(z) is pure imaginary on the intervals (atj, /3j), the numbers Tjk are pure imaginary. 
Furthermore, the matrix r is symmetric and the matrix — it is positive definite (see [9]). 
We now define the Riemann theta function associated with r as 

0(s) = e 2wiim ' s)+7Ti(m ' Tm) , s e C 9 , (4.10) 

m£Z3 

where (m,s) = Y^%X m j s 3- Because the quadratic form i(m,rm) is negative definite, the 
sum in (4.10) is absolutely convergent for all s G C 9 , and thus 8(s) is an entire function in 
C 9 . Notice that the theta function is an even function and satisfies the periodicity properties 

0(s + ej) = 6{s), 6(s + Tj) = e - 2 ™i-™"6{s) (4.11) 

where ej = (0, . . . , 1, . . . , 0) is the j th canonical basis vector in C 9 , and Tj = Ttj. 
Introduce now the vector valued function 

u(z)= [ w, for zeC\(aJ,), (4.12) 



</ 



where u = (ui, . . . , u) g ) is defined in (4.6), and the contour of integration lies in C\(ai, j3 q ) on 
the first sheet of X. Notice that u(z) is well defined as a function with values in C 9 /Z 9 except 
on the interval (ai,{3 q ), where it takes limiting values from the upper and lower half-planes. 
Introduce also the function 

i / _ \ 1/4 

7w=n(f^) (-is) 

with cuts on J°, taking the branch such that 7(2) ~ 1 as z — > 00. It can be seen that, on the 

1 

7 

the function 7 + 

7 

Define the numbers Xj clS 

Xj e (ft, a j+1 ), *y( Xj ) - -J— = 0. (4.14) 

7 \ x j) 

Define the vector of Riemann constants 

<z-i 

K = -^ifij) (4-15) 



first sheet of X, the function 7 — ^ has exactly one zero in each of the intervals ((3j,aj + i), 
and is non-zero elsewhere, and that the function 7 + - has no zeroes on the first sheet of X. 
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and the vector 



Then 



q-l 

d = —K + u(xj). 
j'=i 



(4.16) 



d(u( Xj )-d) = for je{l,...,q-l}, (4.17) 

and {xj}? =1 give all the zeroes of the function 6{u{z) —d). In addition, the function 9(u(z)+d) 
has no zeroes on the first sheet of X. 

Finally, for j — 1, . . . q, introduce the functions 

\ 2/3 




p(t)dt 



^■(*) = - 



3tt 



p(t)dt 



2/3 



1 



(4.18) 



«(oo) + d) 7 (z) + *(«(*) + if + d) 



0(u(oo) + ^ + d) 



6(u(z) + d) 



M 2 {z) 



9{u{oo) + d) 7(2) - 7(^)~ 1 0(u(z) 



2tt 



d) 



(4.19) 



9(u(z) - d) 



0(u(oo) + §f + d) 2 

Notice that .Mi and M. 2 depend quasiperiodically on N, thus are 0(1) as N — > oo. 

The asymptotics of the normalizing constants in equation (1.5) and of the recurrence 
coefficients in equation (1.6) are presented in the following theorem. 

Theorem 4.1. (Asymptotics of recurrence coefficients) Let V(x) be a real analytic func- 
tion satisfying (1-4) which yields a regular equilibrium measure (2.5), and let {P n }^ ^ e ^ e 
system of orthogonal polynomials defined according to (1-5). Then as N — >■ oo ; the normal- 
izing constants in (1.5) and recurrence coefficients in (1.6) admit the following asymptotic 
expansions. 



h N = 



h 



N-l 




oo) + d)6(u(oo) - §f -d) 
G(u(oo) - d)G(u(oo) + ^ + d) 



1 + ° { N 



8Nne 



(w(oo) - d)9(u(oo) 



~~2t 



d) 



9(u(oo) + d)9(u(oo) + §^ - d) 



1 + 



N 



9(u(oo) + d) 2 9(u(oo) - §f - d)9(u(oo) + §f - d) 
9{u{oo) - d) 2 9(u{oo) + §f + d)9{u{oo) - §f + d) 



'N-l 



EJ=iO?-«?) , I V0(^(oo) + *t-d) V9(u(oo) + ^t + d) 
I 0( u(oo ) + 0* - d) " ( M(oo ) + ^ + d ) 



2EJ=i(/3i-«i) 

V0(u(oo) + &) V0(w(oo 



d) 



0(u(oo) + 6) 0(u(oo)-d) 



, u'(oo) 



(4.20) 



(4.21) 



1 

N, 
(4.22) 



(4.23) 
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where V6 is the gradient of 9, 

U'(OO) = (Cl,,_l, C2,,_l, • • • , Cg-lrf-l), (4.24) 

and the numbers Cjk are defined in (4-8). 

Notice that, up to the lattice scaling factor N in the normalizing coefficients, these asymp- 
totics are similar to the results obtained in [6] for continuous orthogonal polynomials. 

The remaining theorems in this section present pointwise asymptotics of the polynomials 
Pn{z) in various regions of the real line and complex plane. 

Theorem 4.2. (Asymptotics of Pn{z) in voids) Let K C C be a compact set on the complex 
plane such that K does not intersect with the support of the equilibrium measure v> . Then 
for any z G K , we have that 

P N ( Z ) = e Na{z) [Mi(z) + 0(N- r )] . (4.25) 

The error term O^N^ 1 ) is uniform in K . 

The function e Ng ^M.i(z) is analytic in a neighborhood of any compact subset of any 
void, thus this formula gives asymptotics of Pn(x) for a; in a void. In particular, notice that 
this function has no zeroes in the exterior intervals (— oo, a±) and (j3 q , oo), and at most one 
zero in any other void. 

Theorem 4.3. (Asymptotics of Pn(z) in bands) Let K be a compact subset of the interior 
of 1°. Then for any point x G K, we have that 

P N (x) = 2e% (v{x)+l) [Re (e iN ^ {x) M 1+ (x)) + O(N^)} , (4.26) 

where Aii+(x) refers to the limiting value of the function Aii(z) from the upper half-plane, 
and 

(f)( x ) ■= / p (t)dt. (4.27) 

J X 

The error term 0(N~ l ) is uniform in K . 

Theorem 4.4. (Asymptotics of Pn(z) in saturated regions) Let K be a compact subset of 
I + . Then the exists e > such that, for any point x G K , we have that 

P N ( X ) = e NL W 2sin(Nirx) (im (e^M 1+ (xf) + C^iV"- 1 )) + 0{e' Ns ) , (4.28) 

where Aii+(x) refers to the limiting value of the function M.\{z) from the upper half-plane, 
and 

L(x):= / log \x - t\p(t)dt. (4.29) 

J ai 

Both of the error terms, O^N^ 1 ) and 0(e~ Ne ), are uniform in K. 

The remaining theorems in this section use the Airy functions Ai and Bi (see, e.g. [13]). 
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Theorem 4.5. (Asymptotics of Pn{z) at band-void edge points) Let j G A v U {q}, so that 
the point (3j is the right endpoint of a band and the left endpoint of a void. Then there exists 
e > such that, for \z — f3j\ < e, 

p N ( z ) = e %{v{z)+i) jjvV6^( z )i/4 Ai(iV 2/3^ (z)) ^ e ±^ M ^ + eT !^ M ^ z) + 0(jV -i)j 

-N- 1 ^p.(z)- 1 ^M , (N 2 ^ p .(z)) e ± ^M 1 {z) - e^^M 2 {z) + O^ 1 ) j 

(4.30) 

for ±lmz > 0. 

Let j G ^U{0}, so that the point etj+i is the left endpoint of a band and the right endpoint 
of a void. There exists e > such that, for \z — cxj + i\ < e, 



P N {z) = e^M+U {iV 1 /^ aj+ij ( 2 )i/4Ai(iV 2 /3 V;Qj+i(2)) 



e 2 



M 2 {z) + 0{N- 1 )^ 



-iV- 1/6 ^ +1 (^)- 1/4 Ai'(iV 2 / 3 ^ 3+1 (^)) e^M^z) + e^^M^z) + 0(N~ X 



for ±Im z > . 



(4.31) 



} 



Theorem 4.6. (Asymptotics of Pn(z) at band- saturated region edge points) Let j G As- 
Then the point f3j is the right endpoint of a band and the left endpoint of a saturated region. 
There exists e > such that, for \z — j3j\ < e, 

i in i N N -, 

■e^M^) + e T ^M 2 (z) + 0(N ) 



P N (z) = e*W'W{N 1 '*il,f ij (z) 1 ''B 1 (z 

- N- 1 / 6 ij f}] ( z y 1 / 4 B 2 (z) \e ± ^ £L Mi{z) + e^^M 2 (z) + O(N^) 



(4.32) 



(4.33) 



for ±Imz > 0, where 

B^z) = cos(N7rz)Ai(N 2/3 ilJ l3j ( Z )) + sm(Nir Z )Bi(N 2/3 ^p.(z)), 

B 2 (z) = cos(Ntt Z )M'(N 2/3 iIj i3] ( Z )) + sm(Nir Z )Bi'(N 2/3 ilJ l3] (z)). 

The point ctj+i is the left endpoint of a band and the right endpoint of a void. There exists 
e > such that, for \z — etj+i\ < e, 



P N {z) = e T(^W+0 {N^^zf^z) [e^M^z) + e^M 2 (z) + OiN- 1 ) 



-N-^.^z)- 1 /^) [e^^-M^z) - e^^M 2 (z) + C^iV" 1 )] ] 

/or ilmz > ; where 

B 3 (z) = cos(N7cz)Ai(N 2 / 3 ^ a]+1 (z)) - sm(Nirz)Bi(N 2 / 3 ^ aj+1 ( Z )), 
B 4 (z) = cos(Nirz)Ai'(N 2 / 3 ip a . +1 (z)) - sm(Nirz)Bi'(N 2 / 3 tp aj+1 (z)). 



(4.34) 
(4.35) 



Remark: Although the above theorems are presented for real analytic potential V(x), these 
results may be extended to potentials which are continuous and piecewise real analytic, 
assuming that the points of non-analyticity lie strictly within saturated regions and voids. In 
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this case the preceding results hold, and the asymptotic solution to the associated Riemann- 
Hilbert Problem does not require local analysis near the points of non-analyticity (see [3]). 



Before continuing with the proofs of these theorems, we would also like to remark that the 
results obtained in this paper match the results obtained in [1] for polynomials orthogonal 
on a lattice which sits inside a finite interval. Consequently, many corollaries discussed in [1] 
also apply to infinite lattices. In particular, the authors of [1] discuss the particle statistics 
of the discrete orthogonal polynomial ensemble in different regions of a finite interval of the 
real line, which are based on asymptotic properties of the associated orthogonal polynomials. 
The results of this paper imply that their results may be extended to discrete orthogonal 
polynomial ensembles on an infinite (regular) lattice. Of particular interest may be the 
discrete sine kernel as the scaling limit of the reproducing kernel in the interior of bands, 
the Airy kernel as the scaling limit of the reproducing kernel at band end-points, the Tracy- 
Widom distribution for the location of the left- and right-most particle, and exponential 
estimates for all correlation functions in voids and saturated regions. 

The rest of the paper is organized as follows. In Section 5, we reformulate the orthogonal 
polynomials (1.4) as the solution to an interpolation problem of complex analysis. In Section 
6, we reduce the interpolation problem to a Riemann-Hilbert Problem which can be solved 
by steepest descent analysis, which is done in Sections 7-12. Finally, in Section 13, we give 
proofs of the preceding theorems. 



5. Interpolation Problem 

We will evaluate the asymptotics of the discrete orthogonal polynomials described above 
via a steepest descent asymptotic analysis of a Riemann-Hilbert problem. To that end, we 
first introduce the following interpolation problem. 

Interpolation Problem. For a given N = 0, 1, . . ., find a 2 x 2 matrix- valued function 
Pn(z) = (PN( z )ij)i<i,j<2 with the following properties: 

(1) Analyticity: Pn(z) is an analytic function of z for z E C \ L N . 

(2) Residues at poles: At each node x G L N , the elements P N (z) u and Pjv(z)2i °f the 
matrix Pjv(^) are analytic functions of z, and the elements Pjv(^)i2 an d Pn{ z )22 
have a simple pole with the residues, 

Res P N (z) j2 = w N (x)P N (x) jl , j = 1, 2. (5.1) 

z=x 

(3) Asymptotics at infinity: There exists a function r(x) > on Ln such that 

lim r{x) = 0, (5.2) 

x— s-oo 

and such that as z — > oo, Pn(z) admits the asymptotic expansion, 



P„(*)~(l + ^ + |i + ...V* N z ° N ), zeC\ 



[j D(x,r(x)) 

.x£Lm 



(5.3) 



where D(x, r(x)) denotes a disk of radius r(x) > centered at x and I is the identity 
matrix. 
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It is not difficult to see (see [4] and [1]) that the IP has a unique solution, which is 

P ( 7 \-( P n( z ) C(w n P n )(z) \ 

N{ ) ~ {(h^^P^z) {h N _ l )^C{w N P N _ 1 ){z)) ' W 

where the Cauchy transformation C is defined by the formula, 

x£L N 

Because of the orthogonality condition, as z — > oo, 

sGLjv x£L n j=0 j=n+2 

which justifies asymptotic expansion (5.3), and have that 

h N = [Pi]i 2 , = [Pxlax. (5.7) 

Furthermore, the recurrence coefficients in equation (1.6) are given by 

7^ = [Pi]i 2 [Pi] 2 i ; fis-i = - [Pi]n- (5.8) 

[f 1J21 

6. Reduction of IP to RHP 

We would like to reduce the Interpolation Problem to a Riemann-Hilbert Problem (RHP). 
Introduce the function 

. sin Nnz 

nw = ^^. (6.1) 

Notice that 

U(x k ) = 0, U'(x k ) = exp (iNnx k ) = (-l) k , for x k = e Lat- (6.2) 
Introduce the upper triangular matrices, 

(1 _«^v(z) ±iNnz\ 
" n(2) i )' (6 ' 3) 



and the lower triangular matrices, 



D± = l-^py^" n(z)J = { n(*)J v -nR^)^ ij • ( 6 - 4 ) 



Define the matrix-valued functions, 



Dl(z) when Imz > 
R« = P N ( Z ) x t + , ~ (6.5) 

' D u z when bnz < 0, K J 



and 



1 D' (A when haz > 

Rjv = P*(*) x ^ + ' " (6.6) 

U_(zJ, when Imz < U. 
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From (5.4) we have that 

1 

W) 



P N (Z) _ ™N(z)P„(z) e ±iN„ Z + C (W N P N )(Z) 



when ± Imz > 0, 



(6.7) 



and 



N[Z) ' - ^^ e^" nW^^P^OWy (6-8) 

when ± Imz > 0. 

Observe that the functions R^(z), R l N (z) are meromorphic on the closed upper and lower 
complex planes and they are two-valued on the real axis. Their possible poles are located 
on the lattice L^. An important result is that, in fact, due to some cancellations they do 
not have any poles at all. We have the following proposition. 

Proposition 6.1. The matrix-valued functions RJ^(z) and R l N (z) have no poles and on the 
real line they satisfy the following jump conditions at x e K: 

R u N+ (x) = R u N .(x)f R (x), f R (x) = (I - 2N ™»(*) ) , ( 6 .9) 



and 



R l N+ (x) = R l N _(x)j l R (x), j l R (x) = (_Lh J) , (6.10) 

\ w N (x) J 



Proof. It follows from the definition of R%(z) that all possible poles of R%(z) are located 
on the lattice L N . Let us show that the residue of all these poles is equal to zero. Consider 
any xt G L^. The residue of the matrix element RJ^ tl2 (z) at Xk is equal to 

Res R« >12 (*) = _ w N(xk)PN(x k ) { _ 1)k + WN{xk) p N{xk) = Q (6 n) 

z — x k y 1 J 

Similarly we get that 

Res Rjv, 22 (z) = 0, (6.12) 

z=x k 

hence R%(z) has no pole at x k - 

Similarly, the residue of the matrix element R l N11 (z) at Xk is equal to 

Res R l N u (z) = - ^M P ^~ l ) k = o. ( 6 .13) 

z=x k N ' UK ' (-1)* WN {x k ) 

In the same way we obtain that 

Res Rjv2i(^) = 0. (6.14) 

Z=Xk 

In the entry R l N2l (z), the pole of the function C(wnPn){z) at z = x k is cancelled by the 
zero of the function H(z), hence R l N21 (z) has no pole at x k - Similarly, R l N22 (z) has no pole 
at Xk as well, hence R l N (z) has no pole at x k - 
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A 



£2' 



£2 



< I 
> 



7^ 

-> — 



CI 



Figure 2. The contour E arising from the hypothetical equilibrium measure 
in Figure 1, dividing an e-neighborhood of the real line into the regions Q± 



and fig. 



Let us evaluate the jump matrices for x G 

r R {x) = BKxy'BKx) 

which proves (6.9). Similarly, 



From (6.5) we have that 



^2isin7V7nr 
n(x) 

1 



-2Ntiiwm{x) 
1 



'+ 



x 



1 

2i sin Nnx 1 



1 

w N (x) 



which proves (6.10). 



(6.15) 

(6.16) 
□ 



To reduce the Interpolation Problem to a Riemann-Hilbert Problem, we follow the work 
[1] with some modifications. Consider the oriented contour E on the complex plane depicted 
in Fig. 2, in which the horizontal lines are Imz = e,0, —e, where e > is a small positive 
constant which will be determined later, and the vertical segments pass through the endpoints 
of saturated intervals. Consider the regions 



nl = {/°ur} x (o,±ie) 

ft* = I + x (0, ±ie) 
bounded by the contour E. Define 



(6.17) 



K^R^K" 1 , 
KjvR^K^ 1 , 



for z G ft±, 
for z G ft±, 
otherwise. 



(6.18) 



where K f 



_(1 * 

If A C C is a set on the complex plane and b G C then, as usual, we denote 

A + b={z = a + b, a G A}. 



(6.19) 



Proposition 6.2. The matrix-valued function Hn(z) has the following jumps on the contour 



R N +(z) = ^N-(z)j R (z), 



(6.20) 
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where 



fl w N (z) 
1 



when z e I U 1° 



, 1 u^QQe*'"' , 

KjvD^K- 1 = [ J aw* n(z) ] 2 e {/" U J } ± ie, 



K JV D' ± (z)K- 1 



n(^)- 1 o 



\-l„±iNnz -rOAT^A„±iNirz J 



K N B l ± {z) 'D^K^ 1 _ \ r 2Nniw N (z)-'e ±tN ™ ^2Nirie= 

when z e (0, ±ie) + /3j or z£ (0, ±ie) + ctj+i for j e A,. 

(6.21) 



7. First transformation of the RHP 
Define the matrix function Tj^(z) as follows from the equation, 

Rjv(z) = e^T^e^H)^ ( 7-1 ) 

where / is the Lagrange multiplier, the function g(z) is described in Section 2, and 03 = 
(0 ^l) * S ^ e ^ird P au ^ matrix. Then T N (z) satisfies the following Riemann-Hilbert 
Problem: 

(1) Tn(z) is analytic in C \ E. 

(2) T N+ {z) = T N _(z)j T {z) for z£E, where 

e N(g_( Z )-±)a 3 j,j e -N(g + ( Z )-±)<j 3 for z( , R 

Jt(z) = I , (7.2) 

\N( g ( Z )-L)* 3jR ( z)e -N {g ( Z )-l)* 3 for 2GS \ M . 



(3) As z — > 00, 



Tjv (z)~I+^ + ^f + .... (7.3) 

z z z 
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Let's take a closer look at the behavior of the jump matrix jr described in (7.2) on the 
horizontal segments of E. We have that 



3t{z) 



when z e I°UI' 



e -NG(z) WN ^ e N(g + ( Z ) +9 -(z)-ir 
e *G{z) 

e~ NG ^ . 

-{2N-K) 2 e- N{9 + {z)+9 - {z) - v{z) ' l) p ng{z) I 11 Z E 

1 ±- 




1 ±- 




i^Y^^j when z = x±ie E {r ± ie} 

e ±NG(z) \ 

_ eT 2W^ e e2iV. \ when ^ = x ± i£ G |j0 ± i£ J 

when 2 = rr ± ie G {/ + ± ie}. 



n(^)- 1 



(7.4) 



8. Second transformation of the RHP 

The second transformation is based on two observations. The first is the well known 
"opening of the lenses" in a neighborhood of the unconstrained support of the equilibrium 
measure. Namely, notice that, for x G 1°, the jump matrix jr(x) factorizes as 

/ e -NG(x) 1 \ _ / 1 OUO lU 1 0\ 

Mx)-y Q e NG{x) j - y e NG(x) ijy-i ojy e -NG( X) x j (g l) 

= j-(x)j M j+(x), 

which allows us to reduce the jump matrix j T to the one 2m plus asymptotically small jumps 
on the lens boundaries. The second observation consists of two facts. Firstly, the jumps 
on the segments I + ± ie behave, for large N , as ±e ±lNnz . Secondly, note that, for x G I + , 
G(x) is a linear function with slope —27ri. With these facts in mind, we make the second 
transformation of the RHP. Let 

' T^zjj+iz)- 1 for zel° x (0,ie) 
T N (z)j-(z) for z G 1° x (0, -ie) 
< T N (z)A + (z) for zel + x(0,ie) 

T N (z)A_(z) for zG/ + x (0, -ie) 
k (2) otherwise. 

f-2lh e ~ iN7VZ ^ ^ - (^^^ ^ 

V "o -2NTtie iN ™) and ~ ^ ™0 2JVfie-* z j' 

(8.2) 

This function satisfies a similar RHP to T, but jumps now occur on a new contour, Eg, 
which is obtained from E by adding the segments (aii — ie, cci + ie), (/3 g — ie, /3 g + ie),(aj +1 — 
ie, + ie), (f3j — ie, f3j + ie) for j G .Ay, see Fig. 3. 



S N (z) = 



where A + (z) = 
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A 



-5»- 



-->- 



/\ 



V 



/\ V 



/\ 



V 



/\ 



/\ \/ 



V 



/\ 



v 



A 



Figure 3. The contour S 5 arising from the hypothetical equilibrium measure 
shown in Figure 1. 



On horizontal segments, we have that 



1 
-1 

1 + 0(e~ 2eN7T ) 0{e N ^~ 2£ ^) 

_ e -NG(z) x 

1 + 0(e- 2eN *) 0(e N ^ G ^- 2E ^) 

e NG(z) 1 

( 1 + 0(e- 2eN7T ) 
\2iNne- N ^)-i-vM) x + 0{e~ 2eN -) 

_ e - iN n j o 

_ e -N{g+(z)+g-(z)-l-V(z)) _ ( ,iNn j 
1 e 7V(2 9 ( Z )-«-V( 2 )) ^ e -2 e JV7r 

1 



-iNfy p N(g+(z)+g-(z)-l-V(z)) 







e iN£lj 



for z £ 1° 
for z — ie £ 1° 
for z + is G 1° 
for z G {/ + ± i£} 
for ze(0j,a j+1 ), jeA s 
for z = x ± G {I~ ± i£} 
for z G (^-,Q!j + i), j G A- 



(8.3) 

By formula (3.5) for the G-function and the upper constraint on the density p, we obtain 
that, for sufficiently small e > and x G (aij, f3j), 



< ^Re G(x ± is) = 2<k P {x) + 0{e 2 ) < 2-ne + 0{e 2 ). 



(8.4) 



This, combined with property (3.7) of the (^-function, implies that all jumps on horizontal 
segments are exponentially close to the identity matrix, provided that they are bounded 
away from the segment (ai,(3 q ). 



9. Model RHP 

The model RHP appears when we drop in the jump matrix js(z) the terms that vanish 
as N — > oo: 



(1) M(z) is analytic in C \ [a±, (3, 
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(2) M+(z) = M-(z)j M (z) for z G [a l7 /3 q ], where 

Jm(z) = { V" 1 °y (9.1) 

-Ms,*** f or z e (^,a j+1 ), 

(3) As z — > oo, 

M(z)~I + ^ + ^ 2 - + .... (9.2) 

This model problem was first solved in the general multi-cut case in [6], and the solution 
is given as follows, using the notation introduced in section 4. 

/ -y( Z ) +1 -^ Z )e(u(z)+^-+d) 7W _ 7 -l( 2 )0(u(^)-gff-d) ' 
M(Z) = F(OO)- 1 2 ' ' 0(u(z)+d) ■ -2i ' e(u{z)~d) \ /g 

\ 2i 9(u(z)-d) 2 e(u(z)+d) 



where 



9(«(oo)+^+d) Q 



F(oo)= l »(«(-,«) mx ;_^ | . (9.4) 

U 0(u(oo)+d) 

The asymptotics at infinity are given as 



M(z) =I+ — + 0(z- 2 ). (9.5) 

z 



10. Parametrix at band- void EDGE POINTS 

We now consider small disks D(aj, e) for j - 1 G A U {0}, and D(f3j, e) for j e AyU {q}, 
centered at the endpoints of bands which are adjacent to a void. Denote 

D=( |J D(a 3 ,e)\\jl \J D(f3 j: e)\. (10.1) 
y-ieA,u{o} / \jeA v u{q} J 

We will seek a local parametrix Un(z) defined on D such that 
(1) 

XJn(z) is analytic onD\ Eg. (10.2) 

(2) 

U JV+ (z)=Uj V _(z)j s (z) for zeDDZs- (10.3) 

(3) 

Ujv(^) = M(z) (I + OiN- 1 )) uniformly for z G 9D. (10.4) 
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We first construct the parametrix near j3j for j G A v . The jumps js are given by 

r 



js(z) 



1 for 2 e (ft - £ > ft) 

_ e -NG(z) for * e (ft, ft + ^) 

e 7VG(.) f) f or * G (ft, ft -ie) 

for z G (/9>-,/9j- +e). 



10.5) 



e -JVG( Z ) e JV(ff+(z)+ fl _(a)-V(«)-J) 

e NG ^ 



If we let 



then the jump conditions on become 

Qjv+(z) = Q N -(z)j Q (z) 



10.6) 
10.7) 



where 



W) 



1 
-1 



1 

-1 1 



(10.8) 



for z G (Pj — e, f3j) 
for z G (j3j, (3j + ie) 
\ f\ for z G (3j — ie) 

(j for ze(pj,pj + e). 

where orientation is from left to right on horizontal contours, and down to up on vertical 
contours, according to Figure 3. 

Q N can be constructed using Airy functions. The Airy function solves the differential 
equation y" = zy, and has the following asymptotics at infinity (see, e.g. [13]): 



AK^^e-I-(l-l,- /2 + 0( ^ 



Ai'(z) = 



2^ 



1/4 -2^3/2 

z 1 e 3 



10.9) 



as z — > oo with arg z G (— n + e, n — e) for any e > 0. If we let 

y (z) = Ai(z), yi{z) = uM(uz), y 2 (z) = u 2 k\(u 2 z) 
where u = e^ 1 , then the functions yo,Vi, and y 2 satisfy the relation 

Vo(z) +yi(z) +y 2 (z) = 0. 



;io.io) 



10.11) 
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If we take 



$1 3g) for - 6 (°- 



-yi(^) -1/2(2) 



7Y 

2 

for argz G ( |,7T 



-1/2 

-2/2(3) 2/i 

»« for arg ,e(4o), 



(10.12) 



y (*) 

then $„, satisfies jump conditions similar to (10.8), but for jumps on rays emanating from 
the origin rather than from We thus need to map the disk D(f3j,e) onto some convex 
neighborhood of the origin in order to take advantage of the function $ ri) . Our mapping 
should match the asymptotics of the Airy function in order to have the matching property 
(10.4). 

To this end, notice that, by (2.23), for t G [a,-, j3j], as t — > /3j, 

p(t) = C(f3 J -t) 1/2 + 0((f3-t) 3/2 ), C>0. (10.13) 

It follows that, for x G [(Xj, f3j] as x — > j3j, 
-Pi 

p(t)dt = C(pj - xf /2 + 0((p j - xf' 2 ) C 



Thus, 



3tt 
T 



p{t)dt \ 



2/3 



(10.14) 



;io.i5) 



is analytic at f3j, thus extends to a conformal map from D(/3j,e) (for small enough e) onto 
a convex neighborhood of the origin. Furthermore, 



1>M) = o ; 1/M)>o, 



(10.16) 



thus ^g. is real negative on (j3j —e, j3j), and real positive on (f3j, f3j+e). Also, we can slightly 
deform the vertical pieces of the contour S5 close to j3j, so that 



We now set 
so that 
where 



n S s } = (-e,e) U (-ie,ie). 
Q N (z) =E^(z)$ rv (N 2 ^(z)) 
U N (z) = E^(z)^ rv (N 2 / 3 ^(z))e- N ^ 



V(z) I 



)<T3 



E#(*) = M(z)e ± ^ CT3 L^(^) _1 for ±lmz>0, 



f N-y^\z) 



2v^F 



# 1/8 <(* 



1 i 
-1 i 



(10.17) 
(10.18) 
(10.19) 

(10.20) 
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and we take the principal branch of ip 1 /^, which is positive on (Pj, (3j + e) and has a cut on 
(Pj — e, Pj). The function $ ru (iV 2 / 3 ^(2;)) has the jumps js, and we claim that the prefactor 
is analytic in D(/3j,e), thus does not change these jumps. This can be seen, as 



M + (z)e^ as = M_(z)e 



lll j,N 

* e 2 



0-3 



l 3,N 



T3 



(10.21) 



thus the jump for the function M.(z)e ~ is 



■'j,N 



■'},N , 



or equivalently, 



1 
-1 



e 2 



CT3 



for z G — e, Pj) 



(10.22) 



u 'j,N 

e 2 



(T3 — /VC3 



al j,N 

e 2 



CM 



for ze(pj,pj + e) 



•'j,N , 



''j.iV 



e 2 ^jMe 2 CT3 



for ze(pj-e,pj) 
(l for zG^ftls), 



;i0.23) 



B ■ B • 3iN Q - 

which is exactly the same as the jump conditions for L#. ThusE^(z) = M(z)e ± 2 a3 L% (z)- 1 
has no jumps in D(Pj, e). The only other possible singularity for is at Pj, and this singu- 
larity is at most a fourth root singularity, thus removable. Thus, E^ is analytic in D(Pj,e), 
and Qjv has the prescribed jumps. We are left only to prove the matching condition (10.4). 
Using (10.9), one can check that, for z in each of the sectors of analyticity, ^^(N^ip^ (z)) 
satisfies the following asymptotics as iV — > 00: 



<S>rv(N 2 ^(z)) 



2v^ 



N-e^.(z)- 



1 A . 4>f>A z )~ 

-1 i) A8N 



3/2 



-5 hi 
-7 -1% 



0(N~ 



x e -|^.( 2 )3/2 CT3 



(10.24) 

where we always take the principal branch of -^.(z) 3 / 2 . As such, -^.(z) 3 / 2 is two-valued for 

z G (Pj — e, Pj), so that 



,3/2 



ft 



p(t)dt. 



Notice that, by (3.9), 

/•ft /-ft 
2g±(x)-V(x) = l±2m p(t)dt = l±2m p(t)dt±iQj 

J X J X 

This implies that for x G (Pj — e, Pj), 

[2g + (Pj) - V(Pj)\ - [2g + (x) - V(x)\ = -2m p(t)dt , 

[2g_(Pj) - V(pj)\ - [2g_(x) - V(x)\ = 2m 



(10.25) 



(10.26) 



/■ft 

/ p(t)dt. 

J X 



(10.27) 
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Combining these equations with (10.25) gives 
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{2g ± {fi j )-V{fi j ))-{2g ± (x)-V(x)) 



(10.28) 



This equation can be extended into the upper and lower planes, respectively, giving 



(2<7±(&) - m)) - W*) - V(z)) 
Since, by (10.26), 2g±(/3 j ) - V(0j) = I ± iflj, we get that 



for ± Im z > 0. 



10.29) 



2 , / \3/2 / x ^l- 2 ) ^ , ^7 



2 2 2 

for ±Imz > 0. Plugging (10.24) and (10.30) into (10.19), we get 



(10.30) 



U N (z) = M(z)e ± ^L^(z)- 1 -^N-^^(z)-^ 



1 % 

-1 i 



+ 



-3/2 



48iV 



-5 5i 
-7 -1% 



+ 0(N~ 2 ) 



M(z) 



e N(g(z)-^-± T ^)* ae -N(g(z)-^l-i)a, 



1 + 48iV \Gte^ -1 ) +G{N 



(10.31) 

for ilmz > 0. Thus we have that Ujy satisfies conditions (10.2), (10.3), and (10.4). 
A similar construction gives the parametrix at the aj for j — 1 G A v . Namely, if we let 



4> aj (z) = - 



3tt 

T 



2/3 



p(t)dt 



10.32) 



then ip aj is analytic throughout D(aij,e), real valued on the real line, and has negative 
derivative at ctj. Close to <x,-, the jumps jq become 



3q{z) = < 



f o r z G (a, - e, aj) 

( 1 0\ , , . . 

I y y\ for z G (ctj, (Xj + zej 

(l l) f° r 2 ^ a i — 

for z G (aj, ctj + e), 



(10.33) 



1 
-1 

where orientation is taken left to right on horizontal contours, and up to down on vertical 
contours according to Figure 3. After the change of variables ip a . (and a slight deformation 
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of vertical contours), these jumps become the following jumps close to the origin: 

for ij} a . ,{z) G (-e, 0) 

'1 N 



1 
-1 



1 1 

1 
-1 1 

1 1 
1 



for i) aj {z) G (0,ie) 

for ^.(2) e (0,-ie) 
for ^q 3 (^) e (0,e), 



10.34) 



where orientation is taken right to left on horizontal contours, and down to up on vertical 
contours. These jump conditions are satisfied by the function 



1 
-1 



Then we can take 

\J N (z) = E%(z)$ lv (N 2 ^(z))e 
for z G D(aij,e), where 



■N(g(z)-^1-L)a 3 



(10.35) 



(10.36) 



E a J(z) = M( Z y 



'j-l,N 



^(z)' 1 for ±lmz>0, 



i I N-y^-^iz) 





2^ 



N^ l J 3 \z) 



1 -i 
-1 -% 



10.37) 



is an analytic prefactor. Similar to (10.24), we have that in each sector of analyticity, 

®iv(N 2 / 3 ip aj (z)) satisfies 



1 -% 



-3/2 



48iV 



-5 -hi 
-7 7i 



+ 0(N~ 



^ lv (N 2 ^ aj (z)) = —N-i°>1, aj {z)-l" [ ^ J. J + 

(10.38) 

Once again, we have that, for x G (aj,aj + e), ip aj (x) 3 ^ 2 takes limiting values from above 
and below, so that 

2 x3 / 2 l = +«ri r p(t)dt. (10.39) 

J Ctj 



;i>ai (x) 



± 



In analogue to (10.28), we have 
2 



{2g ± ( aj ) - V(a 3 )) - (2g{z) - V{z)) 



for ±Imz > 0. 



Since, by (10.26), 2g±(otj) — V(aj) — I ± iri, we get that 
-ib,. ( - V" - - - -ai -. ; + 



i) a \zf t2 = -g(z) + + - ± for ±lmz>0. 



(10.40) 



(10.41) 
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Plugging (10.41) into (10.36) and (10.38) gives, as N ->• oo, 



XJ N (z) = Miz^izr^N-e^iz)-^ 



+ 0(N- 2 
= M(z) 



,N( g (z)-^-l^-^)a 3 -N( g (z)-VM-l)* 3 



i -A V^j (f) 

-1 -i) 48iV 



3/2 



-5 -5i 
-7 7i 



T | ^(^)" 3/2 
A8N 



j-l.N , 



-6ie — a — CT 3 



'j-l,iV 



^3 



+ 0(AT" 2 ) 



10.42) 



11. PARAMETRIX AT THE BAND-SATURATED REGION END POINTS 

We now consider small disks D(aj, e) for j — 1 G *A S , and D(/3j, e) for j G *4. s , centered at 
the endpoints of bands which are adjacent to a saturated region. Denote 



d= ( U ^K^))u( U^' £ )j 

Vj-ieA / Vie A / 



We will seek a local parametrix U N (z) defined on D such that 
(1) 

Uat(z) is analytic on D \ Eg. 

(2) 

U w+ (z)=U w „(2) Js ( 2 ) for zGi)n Eg. 

(3) 

Ujv(^) = M(z) (/ + 0(N~ 1 )) uniformly for z G 3D. 
We first construct the parametrix near (3j for j E A s - Let 

Uat(^) = Q N (z)e TlN ^ 3 e- N ^ z) - X ^ 1 -^ i for ±lmz>0. 
Then the jumps for Qtv are 



1 
-1 

-1 
-1 -1 

1 -1 

1 

1 1 
1 



for z G (Pj — s, f3j) 

for z G (0j,pj + e) 

for z G (fa, Pj + is) 

for z G (/3j, /3j — is) 



(11.1) 



(11.2) 
(11.3) 
(11.4) 

(11.5) 



11.6) 
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where orientation is taken from left to right on horizontal contours, and down to up on 
vertical contours according to Figure 3. We now take 



<f> rs (z) 



for arg z G ^— , 7rj 
for arg z G (-7T, 
?W) for arg,G(-|,0). 



11.7) 



V2(z) yi(z) 

y'i{z) yl(z) 

yi(z) -y2(z) 
y[(z) -y' 2 (z) 

yi(z) y (z) 
y'i(z) y' 

Then <& rs (z) solves a RHP similar to that of Qjv, but for jumps emanating from the origin 
rather than from 
Once again, 

2/3 



Mz) = - 



3tt 

T 



(1 - p(t)) dt 



(11.8) 



extends to a conformal map from D(f3j,e) onto a convex neighborhood of the origin, with 

1>M) = ; i^(&)>0, (11.9) 

Again, we can slightly deform the vertical pieces of the contour £ 5 close to /3j, so that 

(11.10) 



We now take 



where 



^.{%e)nE s } = (-e,e)U(-ie,ie) 
Q w (^ = E^(^$ rs (iV 2 / 3 ^ /3j W) 



/ x x '"j,N 

M(:)< : - 



*-°»I#(z) 



1 for ±Iiri2>0, 



2v^ 



iWV^ 4 (,), 



1 i 
1 -i 



(11.11) 



(11.12) 



and we take the principal branch of ^ • The function ^ rs (N 2 ^ 3 ijj^. (z)) has the jumps j 5 . 

Similar to the prefactor at band- void end-points, the prefactor is analytic in D(/3j,e), 
thus does not change these jumps. 

We now check that XJjy satisfies the matching condition (11.4). The large N asymptotics 
of ^ rs (N 2 ^il)p j (z)) are given in the different regions of analyticity as follows: 



<S> rs (N 2 ^ Pj {z)) 



1 



2v^ 

+ 0(N~ 2 ) 



N-e^ p .(zy 



7<T3 



± 



-i 1 ) A8N 



for ±Im2 > 0, 



-hi 5 
1% 7 



(11.13) 
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where we always take the principal branch of ^.(z) 3 / 2 . As such, ^^.{zfl 2 is two-valued for 

x E (f3j — e, /3j), so that 



3/2 



± 



(1 - p(t)) dt = Tni(f3j -x)±m / p(t)dt. 

J X J X 

From (3.9) we have that 

2g±(x) - V(x) = l±2m / p(t)dt = l±2ni / p(t)dt ± ify =f 2mpj 

J X J X 

for x E (f3j — e, f3j). These equations imply that 

(2g ± (x) - V(x)) - (2g±{Pi) - Vfa)) = ±2ni f' p(t)dt. 

J X 

We can therefore write (11.14) as 



(11.14) 



(11.15) 



(11.16) 



± 



(2g ± (x) - V(x)) - (2g ± (P 3 ) - V(f3j)) 



(11.17) 



We can extend these equations into the upper and lower half-plane, respectively, obtaining 



(2g(z)-V(z))-(2g ± (/3 j )-V(/3 j )) 



for ±1mz > 0. (11.18) 



Using (11.15) at x — /?,-, we can write 



2 -^(zf 2 = g{z) -^- l -±^z T 1{ ^ N * ] for ±1mz>0. (11.19) 
Plugging (11.13) and (11.19) into (11.11) gives 



V N (z) = M(z)e^^(z 



X 



± 



-% -1 

-i 1 



± 



2^ 



N-e^.(z) 



-3/2 



48N 



-hi 5 



+ 0(N~ 2 ) 



x e 



M(z) 



2 2 '"°e _ 

-3/2 



7i 7 
1 -6^^ 



48iV 1-6^^ 



+ 0(iV- 2 ) 



We can make a similar construction near cc,- for j - 1 G ^4 S . Let 



V'a, (*) = - 



3tt 

y 



2/3 



(l-p(t))d* 



for ± Im (z) > 0. 

(11.20) 

(11.21) 
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This function is analytic in D(etj,e) and has negative derivative at aj, thus Imz and 
lmtp aj (z) have opposite signs for z G D(aj,e). Then the jumps for Qtv are 



3Q 



1 
-1 

-1 

-1 -1 
1 -1 

1 

1 1 

1 



for z G (aj, otj + e) 
for z G (aj — e, aj) 
for z G (aj, aj + ie) 
for z G (aj, aj — ie), 



11.22) 



where the contour is oriented from left to right on horizontal segments and up to down on 
vertical segments according to Figure 3. After a slight deformation of the vertical contours 
and the change of variables ip a ., these jumps become the following jumps close to the origin: 



Jq(^a,(z)) 



1 
-1 

-1 
-1 -1 

1 -1 

1 

1 1 
1 



for i; a .(z) e (-e,0) 
for ifj aj (z) G (0,e) 
for ijj aj (z) G (-ie,0) 
for i) a (z)e(Q,ie), 



11.23) 



where the contour is oriented from right to left on horizontal segments and down to up on 
vertical segments. These jump conditions are satisfied by the function 



1 
-1 



11.24) 



<I> ls (z) = § rs (z) 
Then we can take for z G D(aj,e), 

U N (z) = M(2)e^^ CT3 L^(2)- 1 $ is (A^ 2 / 3 ^ aj (z))e Tl ^ 2CT3 e- w(9(z) -^-^ ),T3 for ±lmz>0, 
where 



'N'^~l lA (z) 







2v^ 







N^J^z) 



-1 i 
-1 -i 



(11.25) 
11.26) 



We once again have 
$>is(N 2 ^ aj (z)) 



2v^ 
+ 0(N~ 2 



N-e^ a .(zY 



± 



-i 1\ t/j aj (z) 
-% -1 48N 



-3/2 



-hi - 
7i - 



e l N ^) z/2 °* f or ± lmip a .(z) > (so ^ImoO), 

(11.27) 
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and for z G D(aij,e), 

2 , o /2 . . , . , . Viz) I i(Q jN — n) 

Combining (11.25), (11.27), and (11.28) gives 

XJ N (z) = M{z)e^^^{z)- 1 —=N--^ aj {z)--^ 



for ± Im z > 0. 



x 



i 1 / 48iV 



-3/2 



5i 5 
-7i 7 



+ 0(iV- 2 ) 



x e ±iNnza 3e N(g(z)-^~^)a [ie T l -^a 3e ±lf a 3 ^iNnzas e ~N(g(z)-^ - ±)a 3 



= M(z)e~ 



(T3 



1 + 



-3/2 



48iV 



■1 Qi 



M(z 



-1 



48 N \6ie Ttn i> N 



6, rj +o ^) j 

-0(iV- 2 



e + 2 



0-3 



1 



;n.28) 



for ± Im 2; > 0. 

(11.29) 

12. The third and final transformation of the RHP 

We now consider the contour Ex, which consists of the circles dD(aij,e), and dD(j3j,e), 
for j = l,...q, all oriented counterclockwise, together with the parts of \ (IJ[ a j)/^j]) 

3 

which lie outside of the disks D(a,e), D(a',e), D(/3',e), and D(/3,e), see Fig. 4. 



WU ()-() (K) 





Figure 4. The contour Ex arising from the hypothetical equilibrium measure 
shown in Figure 1. 

We let 

{Sat(2)M(2)~ 1 for z outside the disks D(otj, e), D(j3j,e) 
Sat^Uat^) -1 for z inside the disks D(aj, e), D(/3j,e). 

Then ~K N (z) solves the following RHP: 

(1) Xjv(z) is analytic on C \ Ex- 

(2) Xjv(z) has the jump properties 

Xjv+(i) = X N -(z)j x (z) 



where 



M(z)Ujv(z) 1 for z on the circles 
JX{ ' ] ~ * M(z)j s M(z) ; otherwise. 



12.11 



(12.2) 
(12.3) 
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(3) As z — > oo, 

X N (z) + ^ + (12.4) 

z z z 

Additionally, we have that jx(z) is uniformly close to the identity in the following sense: 

{I + O^N^ 1 ) uniformly on the circles 
I + (e" c ^) on the rest of S x , ^ 
where C(z) is a positive, continuous function satisfying (1.4). If we set 

&(z)=j x (z)-I, (12.6) 

then (12.5) becomes 

(0(N~ 1 ) uniformly on the circles 
0(e-°M") on the rest of E x . ^ 
The solution to the RHP for is based on the following lemma: 
Lemma 12.1. Suppose v(z) is a function on T, x solving the equation 

v(u)j x (u) 



;i2.8) 



2m J^ x z--u 

where Z- means the value of the integral on the minus side of T, x . Then 

X N (z)=I-—[ v{u)3 ° M du forze£\Y> x (12.9) 

solves the RHP for X^r- 

The proof of this lemma is immediate from the jump property of the Cauchy transform. 
By assumption 

X N _(z) = v(z) (12.10) 
and the additive jump of the Cauchy transform gives 

X N+ (z) - X N _(z) = v(z)j° x (z) = X N _(z)f x (z) (12.11) 

thus Xn+(z) = ~X.N-(z)jx(z). Asymptotics at infinity are given by (12.9). 

The solution to equation (12.8) is given by a series of perturbation theory. Namely, the 
solution is 

oo 

v(z) =I + ^2v k (z) (12.12) 

k=l 

where 

v k -i(u)jSc(u) 



Vk(z) = -— / v /JAV '- du ; v (z) = I. (12.13) 

2m J Sx z-u 

This function clearly solves (12.8) provided the series converges, which it does, for sufficiently 
large N. Indeed, by (12.5), 

/ C\ k 1 

|f fc (z)| < yj^J i -|- | | ^ or some cons tant C > , (12.14) 
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thus the series (12.12) is dominated by a convergent geometric series and thus converges 
absolutely. This in turn gives 



X N (z) = I + J^X N>k (z) (12.15) 

where 



k=i 



x^) = --L / "tiWIMi. (12 . 16) 

2m /v„ z — u 



In particular, this implies that 

x "~ I + Q ( jv(|*| + i) ) as N ^°° (12 ' 17) 

uniformly for z G C \ Ex • 

13. Proof of theorems 4.1-4.6 

The transformations (6.18), (7.1), (8.1), (12.1) give that, for z bounded away from the 
real line, 

Pjv(z) = K- 1 e^ a3 X N (z)M(z)e N ^-^ a3 K N , (13.1) 
and for z close to the real line but bounded away from the support of the equilibrium measure, 

P N (z) = Kj f 1 e I: T a3 X N (z)M(z)e N{9{z) -^ a3 K N T>l(z)- 1 for ±lmz>0. (13.2) 
Expanding (13.1) or (13.2), we get that 

P N (z) = [Pjv(z)]h = e^W ([M] n [X] n + [M] 21 [X] 12 ) (13.3) 

which, along with (12.17), proves Theorem 4.2. 

The proof of Theorem 4.1 requires only the formulae (5.3), (5.7), and (5.8), and a straight- 
forward large- ,2 expansion of equation (13.1). 

Similar to (13.1), we have that, for any interval J which is contained in and bounded away 
from the endpoints of a band, in some neighborhood of J, we have 



) K N 1 e^ a3 X N (z)M(z)j + (z)e N(9{z) -^ )a3 K N -Dl(z)- 1 for Imz > 

K N 1 e^ a3 X N (z)M(z)jZ\z)e N{ ^ z) '^K N B u _(z)- 1 for lmz<0, ^ 3 '^ 

Expanding the left side of this equation for Imz > 0, utilizing (3.9), and taking limits as z 
approaches the real line, we get that 

P N ( X ) = [P N (x)] n = eTOW) ( e ^^[M 11 ] + (x) + e"^ (a;) [M 12 ] + (x) + OiN' 1 )) 

(13.5) 

where <f>(x) is as defined in (4.27), and the + subscript indicates the limiting value from the 
upper half plane. Notice that [M 12 ]+ = [M n ]_ in this region, and that M_ n (z) = Mn(z). 
This implies that [Mi 2 ] + (i) = [Mn] + (i), and thus we can write (13.5) as 

P N {x) = e^ v(x)+l) (e iN ^ x) [M n ] + (x) + e iN **(*)[M n ] + (x) + 0(N^ , (13.6) 

which proves Theorem 4.3. 
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= e NL{x) 



For any interval J which is contained in and bounded away from the endpoints of a 
saturated region, in some neighborhood of J, we have 

P N (z)^K- 1 e^ as X N (z)M(z)Al 1 (z)e N ^-^ a3 K N B l ± (z)- 1 for ±lmz>0. (13.7) 
Notice that in this region, we can write 

g±(x) = L(x)± l ^-Ti™ (13-8) 

where L(x) is defined in (4.29). Notice also that 2g±{x) — V(x) — I has positive real part. 
Expanding (13.7) for Imz > and taking the limit as z approaches the real line gives 

P N (x) n = e^+W [(1 - e 2 ™ Nx )(M n X n + M 21 X 12 ) + e' N ^+^~ v ^' l \M 12 X n + M 22 X 12 )] 

= e N9 + {x) [(1 - e 2 " JV:E )(M 11 X 11 + 0(N- 1 )) + 0{e~ m )] 

-2i sm(irNx)e^ [M u ] + (x)(l + 0(N~ 1 )) + 0(e- N5 ) \ , 

(13.9) 

which proves Theorem 4.4. 

Similarly, at the turning points Oij and {3j, explicit formulae can be written for in terms 
of explicit transformations in each sector of analyticity of the local parametrix. From these 
formulae and the properties of the ^-function, Theorems 4.5 and 4.6 are almost immediate, 
with Theorem 4.6 also requiring the identities (see, e.g. [13]) 

y 1 (z) = ~(M(z)-im(z)) 

\ (13.10) 

y 2 (z) = --{M(z)+iBi(z)). 
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